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J> , A three-by-three matrix spectral problem for AKNS soliton hierarchy is 

proposed and the corresponding Bargmann symmetry constraint involved in 
Lax pairs and adjoint Lax pairs is discussed. The resulting nonlinearized Lax 
systems possess classical Hamiltonian structures, in which the nonlinearized 
spatial system is intimately related to stationary AKNS flows. These nonlin- 
■ earized Lax systems also lead to a sort of involutive solutions to each AKNS 

soliton equation. 



1 Introduction 



Symmetry constraints have aroused an increasing interest in recent few years due to 
the important roles they play in soliton theory. Such a kind of very successful sym- 
metry constraint method is the nonlinearization technique for Lax pairs of soliton 
hierarchies, including mono-nonlinearization proposed by Cao and Geng || and 
further binary nonlinearization [T7J [TJj]. 

In general, one considers the complicated nonlinear problems to be solved in 
such a way to break nonlinear problems into several linear or smaller ones and then 
to solve these resulting problems. It is following this idea that one has introduced 
the method of Lax pair to study nonlinear soliton equations. The Lax pairs are al- 
ways linear with respect to their eigenfunctions. Nevertheless, the nonlinearization 
technique puts this original object, the Lax pair, into a nonlinear and more compli- 
cated object, the nonlinearized Lax system. It seems to be not reasonable enough, 
but in fact, it provides an effective way, different from the usual one, to solve soliton 
equations. The main reason why the nonlinearization technique takes effect is that 
kind of specific symmetry constraints expressed through the variational derivative 
of the potential. 
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A similar symmetry constraint procedure for bi-Hamiltonian soliton hierarchies 
is presented by Antonowicz and Wojciechowski et al H H [2pJ and bi-Hamiltonian 



structures for the resulting classical integrable systems can also be worked out 
through a Miura map [|]. A connection between these systems and stationary 
flows 0| is also given by Tondo [24] for the case of KdV hierarchy. Because station- 
ary flows may be interpreted as finite dimensional Hamiltonian systems [[J based 
upon the so-called Jacobi-Ostrogradsky coordinates [TO, a natural generalization of 



nonlinearization technique to higher order symmetry constraints is made by Zeng 
[ 27fl [28] for the KdV and Kaup-Newell hierarchies etc. There have also been some 
algebraic geometric tricks, proposed by Flaschka et al [11 1 [|I]] [23], to deal with 
similar nonlinearized Lax pairs called Neumann systems. 

The study of the nonlinearization theory leads to a large class of interesting 
finite dimensional Liouville integrable Hamiltonian systems which are connected 
with soliton hierarchies (for example, see || |15|]). However in the literature, most 
results are presented for the cases of 2 x 2 matrix spectral problems. The present 
paper is devoted to the symmetry constraints in binary nonlinearization for a case 
of 3 x 3 matrix spectral problems. We successfully propose a 3 x 3 matrix spectral 
problem for AKNS soliton hierarchy, motivated by a representation of 3 x 3 matrices 
for the Lie algebra sl(2). Then in Section 3, we consider the Bargmann symmetry 
constraint for the proposed new Lax pairs and adjoint Lax pairs of AKNS soliton 
hierarchy. In Section 4, we analyze the nonlinearized Lax systems, especially the 
nonlinearized temporal systems, and establish a sort of involutive solutions to AKNS 
soliton equations. Finally in Section 5, some remarks are given. 



2 New Lax pairs for AKNS equations 



We introduce a three-by-three matrix spectral problem 



U{u,X) 



where the potential u = (q, r) T . Its adjoint spectral problem reads as 
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Here T means the transposition of the matrix. Our purpose is to generate AKNS 
hierarchy of soliton equations from the above specific spectral problem Q2.1p . To this 
end, we first solve the adjoint representation equation V x = [U, V}. Take 
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and then we have 



[U,V] 



I 2(qc-rb) -2V2(\b + qa) \ 

2 v / 2(ra + Ac) -2 V / 2(A6 + qa) 

V 2V2{ra + \c) -2{qc-rb) J 



Therefore we easily find that the adjoint representation equation V x = [U, V] be- 
comes 

a x = qc — rb, b x = —2Xb — 2qa, c x = 2Xc + 2ra, 
which is equivalent to 

a ix = qci - rbi, b ix = -2b i+1 - 2qa iy c ix = 2c m + 2ra { , i>0. (2.4) 

We fix the initial values 

a = -1, b = c = (2.5) 

and require that 

di\ u =o — h\u=o — c i\u=o = 0, i > 1, (2-6) 

which equivalently select constants of integration to be zero. On the other hand, the 
above equality (2^) gives rise to the recursion relation for determining a,, bi, cf. 

a i+ i = \d~ l (qc ix + rb ix ), 

h+x = -\bix - qai, i > 0. (2.7) 

Q+l ^C-ix TQiii 

This recursion relation uniquely determines infinitely many sets of polynomials 
di,bi,Ci, i > 1, in u,u x ,--- under the requirement ( [2.6|) . The first two sets are 
as follows 

1 11 

ai = 0, h = q, a = r; a 2 = -(gr), b 2 = —j<lx, c 2 = -r x . 

In addition, we have 

CO CO CO 

a 2 + bc= (£ a, A-) 2 + (£ biA-^E aX-*) = 1, 

i=0 i=0 i=0 

because (a 2 + bc) x = ^tr(V 2 ) x = |tr[[7, V 2 ] = and (a 2 + bc)\ u=Q = 1. It follows that 
Q'ii h, Ci, i > 1, are local. 

A direct computation may show that the compatibility conditions of the Lax 
pairs 

<P X = <P tn = V^<P, = V^(u, A) = (X n V) + , n>0, (2.8) 
or the adjoint Lax pairs 

^ = -U T ^ i> tn = -(F W ) T ^, n > 0, (2.9) 
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where the symbol + denotes the choice of non- negative power of A, engenders a 
hierarchy of AKNS soliton equations 

« t =(«) t = A.= (-^')= J L-( ;),„><>, (2.10) 

where the Hamiltonian operator J and the recursion operator L read as 

/0-2\ (Id-rd-'q rd- l r \ 

J -\2 )' { -qd-'q -\d + qd^r J ' 

This AKNS hierarchy is exactly the same as one in Ref. [f[8|| , which also shows 
that the same soliton hierarchy may possess different Lax pairs, even different order 



spectral matrices. Here the operator L* is a hereditary operator | 13|j , and J and JL 
constitute a Hamiltonian pair. 

Finally, we would like to elucidate the other two properties on AKNS hierarchy 
( |2. 10] ). First by Corollary 2.1 of Ref. we can obtain 



V tn = [V^,V], n>0 (2.12) 

when u tn = K n , i.e. U tn — + [U, V^} = 0, n > 0. Second, we can get the 
Hamiltonian structure of AKNS hierarchy 



u tn = K n = J ( ° h n+1 ) = J^, H n = -?—a n+2 , n > 0, (2.13) 
Wi / du n+1 



by applying the trace identity p5| p6| 



3 Binary nonlinearization related to new spectral 
problem 

In order to impose the Bargmann symmetry constraint in binary nonlinearization, 
we first need to compute the variational derivative of the spectral parameter A with 
to the potential u, which is shown in the following Lemma ||12|||]I5 |. 

Lemma 3.1 Let U(u, A) be a matrix of order s depending on u, u x , ■ ■ ■ and a pa- 
rameter A. Suppose that (j) = (<pi, fa, ■ ■ ■ , 4> S ) T , ip — C0i) fa, 4 • • , "0s) T satisfy the 
spectral problem and the adjoint spectral problem 

(j) x = U(u,\)(f>, tp x = -U T {u,\)%p, 

and set the matrix V = (pifj T = {4>kfa)sxs! then we have the following two results: 

(i) the variational derivative of the spectral parameter A with respect to the potential 
u may be expressed by 

8X tr(V§^) 



5u -JZoHVmdx' 
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(ii) the matrix V is a solution to the adjoint representation equation V x = [U,V], 
i.e. V X = [U,V]. 

Following ( |3.1|) , we have the variational derivative of the spectral parameter for 
the spectral problem (|2.1|) and the adjoint spectral problem (|2.2| ) 



(3.2) 



5X y/2 5X a/2 , , i , \ 



where E = 2/^ o (0i'0i - fcip^dx. 

Let us introduce N (N >l) distinct eigenvalues Xj, 1 < j < N, and denote by 

<j>W = (0 lj ,0 2j ,0 3j ) T , 1>® = (iJ; lj ,i> 2j ,ij 3j ) T , 1 < j < N, 

the eigenfunctions of (|2.S| ) and the adjoint eigenfunctions of (|2.9|) , i.e. 

= £/(«, Xj)(fr( j \ ^ = -U T (u, X 3 )^\ l<j<N, (3.3) 

^ = 7( n »(«, X 3 )^\ 1$ = -{V^) T {u, X^\ l<j<N. (3.4) 
Now we make the Bargmann symmetry constraint 



Kn = J 



8u 



N ^ 



5w ' 



(3.5) 



where Ej = 2/^ o (0i i '0y - <j) 3 jip3j)dx, 1 < j < N, and 1 < j < A/", are any 
nonzero constants. By ( |3.2|) , this symmetry constraint becomes 

i, f ^{fcj^j + 03^) ^ 



from which we get the following explicit expression for the potential u 



u 



/(ft, ft, ft; ft, ft, ft) = V2 ( < ™ > + < ™ > ) . (3. 6) 



Here and hereafter, < •, • > denotes the standard inner product of H N and 



B = diagQ/i,- ■ - ,/ijv), 



Pi 
Qi 



Wil, <Pi2, "', <PiN) 



,1 = 1,2,3. (3.7) 



The substitution of ( p.6|) into the spatial system (|3.3|) and the temporal systems 
fl3.4|) for n > yields the nonlinearized spatial system: 





f 0y \ 




02j 
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f 1% \ 











' <pij 



, J" = 1,2, 



(3.8) 



-U (f, Xj) | ^ | , j = 1,2, 

^3j 



5 



and the nonlinearized temporal systems for n > 0: 





( <h \ 
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•(v (n) ) r (/,A 




(3.9) 



1,2,---,N. 



~M is a 



It is obvious that (|3.8|) is a system of ordinary differential equations and 
hierarchy of partial differential equations. 

Suppose that Z is an expression depending on u and its differentials. From now 
on we use Z to denote the expression of Z depending on P i: Q t , 1 < i < 3, and their 
differentials after substituting ( |3TB| ) into Z, and use Or(Z) to denote the expression 
of Z only depending on Pi,Qi, 1 < z < 3, themselves after substituting ( |3.8j ) into 
Z sufficiently many times. Therefore ( |3.9| ) may be transformed into the following 
systems for n > 0: 
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j = l,2,--.,N, 



j = l,2,---,N, 



(3.10) 



which are all ordinary differential equations with an independent variable t n because 
the matrices Or(V^ n \f, A 3 -)), n > 0, 1 < j < N, only depend on P { , Q i5 1 < i < 3. 

We would like to discuss the integrability on the nonlinearized spatial system 
Q3.8| ) and the nonlinearized temporal systems ( |3.10| ) for n > in the Liouville sense 
2lfl . We shall utilize the symplectic structure uj 2 on R a7V 



3 N 3 

Y Y H d( t>ij A dipij = Y( Bdp i) A d Qi 

i=0 j=0 



i=0 



by which one can define the corresponding Poisson bracket for two functions F, G 
defined over the phase space R 6Ar 



{F, G} = u 2 {IdG, IdF) = ou 2 (X G , X F ) 



, i , : <)c,j <)<:>,: ()u,j i)i.' u 
*, dF .dG dF 



dQ l 
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where IdH = Xh denotes the Hamiltonian vector field with energy H determined 
by 

u 2 (X,IdH) = u 2 (X,X H ) = dH(X), X G T(R 6Ar ), 
and the corresponding Hamiltonian system with the Hamiltonian function H 

x = IdH(x) = dx(IdH) = u 2 {IdH, Idx} = {x, H}, xe R 6JV , (3.12) 

which possesses an explicit formulation 

Note that there are some authors who use the other Poisson bracket {F, G} = 
u 2 (X F , X G ). As remarked by Carroll ||, it doesn't matter of course but each type 
has many proponents and hence one must be careful of minus signs in reading various 
sources. The notation we accept here is the Arnold's one Hj. 



Theorem 3.1 The following functions 



i- ) y.o^dj. i<j<N, 



(3.14) 



i=i 



are all integrals of motion for the nonlinearized spatial system (\3.8j ). Moreover they 
are in involution under the Poisson bracket (°J.ll) and independent over the region 



n = {R m \ <p iJ} ^ G R, YMi + $i) ± °> 1 ^ 3 < N}. 



Proof: Let 



We can first find that 



V{\j) = (0jy^y)fc,i=iA3 , 1 < j < N. 



tr(V^)). 



On the other hand, by Lemma [37L] we know that V(\j) satisfies 

V{^j)x = [U{u, Xj),V(Xj)] 

when (|3.8|) holds, and thus 

= tr[£/(S,A i ),V'(A i )] = 0, 

which shows that Fj, 1 < j < N, are all integrals of motion for the nonlinearized 
spatial system (|3.8|) . In addition, it is very easy to prove that 

{F kl F l } = 0, l<k,l<N, 
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which means Fj, 1 < j < N, are in involution. It is also obvious that gradFj, 1 < 
j < N, are everywhere linear independent over Q by observing that 



/ ^ 



( & ) 

\ dr Pu / k,i=i, 
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( e£k ^ 
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i = 1,2,3, 



>t2 



, » = 1,2,3. 



4>iN j 



The proof is completed. I 

Throughout our paper, we assume that 

A = diag(Ai, A 2 , • • • , Aat) 



(3.15) 



If all elements Zij, 1 < i,j < s of a given matrix Z = (Zij) sxs are polynomials in 
A, i. e. Z^ = Yl™=o ^ijk^ k , for convenience of presentation, we define a new matrix 
called Ma(Z) as follows 



sNxsN- 

k=0 

Moreover we often accept compact forms, for example, 

d , d d 



(3.16) 



dR 



Y , {P t , H} = ({0a, H}, ■ ■ ■ , {& N , H}) 1 , i = 1,2,3. 



(3.17) 



Theorem 3.2 We have the explicit integrals of motion for the nonlinearized spatial 
system 

{ F 1 = -8(< P 1 ,BQ 1 >-< P 3 , BQ 3 >), 
F m = 4E^ 1 [(< A i ~ 1 P 1 ,BQ 1 >-< A i ' 1 P 3 ,BQ 3 >)x 
(< A m - l - 1 P 1 ,BQ 1 >-< A m -'- l P^ BQ 3 >) 
+2(< A l - 1 P 1 ,BQ 2 > + < A i ~ 1 P 2 ,BQ 3 >)x 
(< A m - i - 1 P 2 , BQ X > + < A" 1 "*- 1 ^, BQ 2 >)] 
-8(< A m - x P u BQ x >-< A m - l P 3 , BQ 3 >), m > 2, 

where Pi, Qi, B are defined by ( \3. T) . Moreover they constitute an involutive system 
together with Fj, 1 < j < N, under the Poisson bracket ( \3.11\ ), i.e. 

{F k , Fi} = {F m , Fj} = 0, m, k, I > 1, 1 < j < N. 
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Proof: We assume that 



q = V2{< Pi, BQ 2 > + < P 2 , BQ 3 >), f = y/2(< P 2 , BQ 1 > + < P 3 , BQ 2 > 
So = -1, b = c = 0; 

a l+1 =< A'P,, BQi > - < A l P 3 , BQ 3 >, % > 0, 

k+i = \/2(< A'Pi, PQ 2 > + < A J P 2 , PQ 3 >), i > 0, 
Ci+i = v / 2(< A*P 2j PQi > + < A'P 3 , BQ 2 >), i > 0. 



(3.18) 



Further we choose that 



-2A V2q 
T ! v^r V2q \ , V 
v^r 2A 



' Ta \ 

v^c V2b 
K v^c -2a j 



where a, 6 and c are defined by 



i=0 i=0 i=0 

It may be shown that when the nonlinearized spatial system ( |3.8| ) holds, we have 

V x = [U, V], i.e. a x = qc — rb, b x = —2Xb — 2qa, c x = 2Ac + 2fa. 
Therefore we can compute that 

F* ■= (^(V")) x = ±tr(V% = hrp,^] = 0. 
On the other hand, we have 

oo m 

F = 4(a 2 + be) = F m X- m , F = 4, F m = 4]T(a i a m _ i + M™-*), m > 1. 

m=0 i=0 

Hence F m , m > 1, are all integrals of motion for the nonlinearized spatial system 



Now we turn to the involutivity of integrals of motion. We take 

V (n) (A) = (X n V)+ 
and construct a temporal system for n > 



Pi 

P-2 
P 3 



Pi 

M A (X>^) I P 2 
P 3 



Qi 
Q2 



Qi 

-(M A (V*»>)) T | Q 2 
Q3 



(3.19) 



where M A (V {n) ), n > 0, are determined in the way of ( |3.16|) . We can first prove that 
when this system Q3.19 ) holds, we have 
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Therefore F m , m > 1, are also integrals of motion for the system fl3.19| ). Secondly, 
we can verify that 



V B-% J 



p-l dF 
D dQ 2 



9Q3 / 



I B-Hx(VJLV) \ 
B-HriVjL-V) 



I B-hriV^V) \ 

b-Mv^v) 



m=0 



m=0 



( Q1 \ 

Q2 



-m—1 



(Pl\ 
P2 



A 



-m— 1 



V b-Mv^-v) J 

These two equalities show that the system ( |3.19 ) for n > are all Hamiltonian 
systems with Hamiltonian functions —F n+ i. Therefore 

{Fm+i, -F n+1 } = -^-F m+1 = 0, m, n > 0, 

(Ion 

which shows the involutivity of F m , m > 1. In addition, it is easy to get that 

w,*.}=£rt-'(|£-' SFk 



il 



HI 



0, k > 1, I < I < N, 



i=Q ' dvpu ucpu 
noting the particular form of F^, k > 1. The proof is finished. I 



Because we have V 2 (X j ) = tr(y(A j ))V r (A i ), 1 < j < N, V 3 = §tr(F 2 )V ', we 



cannot obtain new integrals of motion of the nonlinearized spatial system ( ^.8[ ) from 
the trace of other power of V(\j) and V. In addition to this, it is interesting to 
observe that the determinants of the matrices V(Xj), 1 < j < N, and V are all zero. 

The nonlinearized spatial system is easily rewritten as an Hamiltonian system 



P lx = {P u H} 



-B Qi 



f)J-f 

{Q l ,H} = B- 1 ™, i 



1,2,3 



1^ + h F l 



with the Hamiltonian function 
H = 2(< AP u BQ l >-< AP 3 , BQ 3 >) 

-2(< P 1 ,BQ 2 > + < P 2 , BQ 3 >)(< P 2 , BQ 1 > + < P 3 , 5Q 2 >) 

Thus it possesses the following 3iV involutive integrals of motion 

Fj, 1 < j < N, F m , 1 < m < 2N. 

In some special cases, they may be shown to be independent at least on certain 
region. 
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Theorem 3.3 When N = 1,2, the determinant of the matrix 

( {gradpF^ (gradp^f (gradp^f \ 



D(N) 



(grad Pl F N ) T (gradp 2 F N ) T 



"Nf {9 r aap 2 fN)~ (grad Pa F N ) rj 
(grad Pl Fi) T (gradp 2 Fi) T (gradp^) 7 



\ (grad Pl F 2N ) (gradp 2 F 2N ) T (grad Ps F 2N ) T J 



with gradp G 



dG 



dG \T 



, 1 < i < 3, is not always zero. Thus the integrals of 



motion Fj, 1 < j < N, F m , 1 < m < 2N, are independent at least on some region. 



Proof: We have 

D(l) = {lGfaMll + 16031^21^31 - 320x^^^21^31 

+32021^11^21^31 - 32031^11^21^31 _ 640 2 l'0ll'03l)A t l 

and 

D{2) = (-256A? + 256A^ - 768A1A2 + 768A?A 2 ))U?/i 2 , 
when we choose 

011 = 1, 012 = 0, 021 = 0, 022 = -1, 031 = 0, 32 = 0, 
Ipll = 0, 1pl2 = -1, 1p21 = 1, ^22 = 0, -031 = 0, ^32 = ~1- 

These are consequences of computation by the computer algebra system MuPAD. 
But they may also be shown by some direct computation. The proof is completed. I 

According to the above theorem, the nonlinearized spatial system ( |3.8j ) is Liou- 
ville integrable on some region of the phase space, when N = 1,2. 



4 Involutive solutions 



The aim of this section is to discuss some properties on the nonlinearized spatial 



system ( |3.8|) and the nonlinearized spatial system ( |3.10|) and to establish a kind of 
involutive solutions with separated variables for AKNS soliton equations. 

Lemma 4.1 When <fi = (0i, 2 , 03) r and if) = (ipi, -02, "03) T satisfy the spectral prob- 



lem (\2.1j ) and the adjoint spectral problem ( ]2.2\ ), we have 

L j 02^1 + 03^2 \ _ A ( 02^1 + 03^2 \ j 
01^2 + 02^3 



02^1 + 4>3^2 
<t>li>2 + 02^3 



r 
Q 



(4.1) 



where I is an integral of motion for ( \2. 1\ ) and ( \2.2\ ). 
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Proof: From the spectral problem ( j2.1p and the adjoint spectral problem (|2.2[ ), we 
can find that 



9(01^1 - h^) = V2q(<fr 2 A + 03^2) - v / 2r(0 1 ^ 2 + 0203)- 

This yields 

9~ 1 [-g(02'01 + 0302) + K0102 + 0203)] = ~ 03*03) + 



where 7 is an integrals of motion for ( |2.1| ) and ( |2. 2| ) . The relation (|4. 1|) follows from 
the above equality. I 

We recall that Z denotes the expression of Z depending on P i: Q i} 1 < % < 3, and 
their differentials after the substitution of ( |3.6| ) and into Z, and that Or(Z) denotes 
the expression of Z only depending on Pi,Qi, 1 < i < 3, themselves after the 
substitution of ( |3.8| ) into Z sufficiently many times. A general result on Sj, bi, Ci, i > 
1, is given in the following theorem. 

Theorem 4.1 VKe have the explicit expressions for a m , b m , c m ,m > 1: 

a m+1 = 5^ /,(< A^P^BQr >-< A^P^BQs >) - I m+1 , m > 0, (4.2) 

4 = 

m 

b m+1 = v / 2E^(< A m ~ l P u BQ 2 > + < A m -*P 2 ,BQ 3 >), m > 0, (4.3) 

i=0 
m 

£ m+1 = v^E J i(< ^-^2, BQi > + < A^Pg, PQ 2 >), m > 0, (4.4) 



i=0 



provided that the nonlinearized spatial system ( \3.Sj ) is satisfied. Here I m , m > 0, are 
defined by 

m 

I = l,I m = Y,d n £ F h ---F in , m>l, (4.5) 

ra=l 

«lH h«n=TO 

where the constants d n , n > 0, are determined recursively by 

r di = -i ^ = its, 

I W — —I v rt_1 W ^7 — 1(7 — I v n ~ 2 ///I n > ^ v • / 

u n 2 <£-'i=l Uj i u "n—i 4 U/ n— 1 g Z^j=l u, i Uj n—i—li 11 _ J j 

and iae functions F m , m > 1, are azwen fry t \3.1% ). 



Proof: By the recursion relation (|2.7|) and Lemma ( |4.1| ), we can obtain that 



c m+1 \ = ~ m ( r \ = ~ ( V2(< P 2 .BQ 1 > + < P 3 .BQ 2 >) 
b m+ i J \q) \V2{<P 1 ,BQ 2 > + <P 2 ,BQ 3 >) 

™ / v^(< A m ~*P 2 , SQx > + < A m-i P 3 , BQ 2 >) \ 
^ 4 ^ v^(< A^Px, PQ 2 > + < A m -P 2 , BQ 3 >) )' m - U ' 
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where Jo = 1 and Jj, 1 < i < m, are integrals of motion for the nonlinearized spatial 
system (|3.8|). Now we compute a m+ i, m > 0, by b mx = —2b m+ i — 2qa m , m > 0. 
Noting (|3.8|) , we have for m > 



m— 1 



b mx = J2^(< Am ^ lp ^^ B Q2> + <A m - 1 - 1 Pi,BQ 2x > 

+ < A m ~ i ~ 1 P 2x , BQ 3 > + < A m - l - l P 2 ,BQz x >) 

m— 1 

= J2 -2^ m_i Pi, Sg 2 > + < A m ~ i ~ 1 P 1 , -V2qBQ 1 > 

+ < V2qA nM P 3 , BQ 3 > + < A m - i - 1 P 2 , -2ABQ 3 >) 

m— 1 

- -26 m+ i -2g(E I i (< A m ~ i ~ 1 P 1 , BQi >-< A m -^ 1 P 3 ,BQ 3 >)-I m ), 

i=0 

from which ( |4.2| ) follows. 

In the following we determine the integrals of motion I m , m > 0, by a relation 



a 2 



+ 6c = 1, a = ^ajA l , b = ^2bi\ \ c = ^CjA \ 



i=0 i=0 «=o 

which gives rise to 

m— 1 

2a m = ^ (a;a m _; + ^c m _i), m > 2. (4.7) 



i=l 

First from Si = we have 



h =< P l5 SQx > - < P 3 , BQ 3 >= -]-F u 



which shows d\ = — |. Now we suppose m > 2. At this moment, we have by ( |4.7| ) 

m— 1 

2 ^ ] li&m—i 2/ m 



i=0 

m— 1 i— 1 m—i—l 

E (X/ Ik&i-k — E hO"m-i-l — Im- 

i=l fc=o Z=0 

m— 1 i— 1 m—i—l 

+ E E E hcm-i-h rn>2, 

i=l fc=0 z=o 



where eij, 6j, Cj, i > 1, are given by fl3.18|) . After interchanging the summing in the 
above equality, i.e. 

m— 1 i— 1 m-2 m— 1 m— 1 m—i—l m— 2 m-l-Z m— 1 i-1 m— i-1 m-2 (m-2)-k m-(Z+l) 

EE = EE>EE=EE>EEE=E E E . 

i=l fc=0 fc=0 i=fc+l i=l Z=0 Z=0 i=l i=l fc=0 Z=0 k=0 1=0 i=k+l 
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we may arrive at 

m-l m—2 (m-2)-fc m-(£+l) 

— 2I m = E Idm-i + E E E Ikhifli-k&m-i-l + bi-kC m -i-l) 

i=l fc=0 i=0 i=fc+l 

m—2 m— 1 m-2m-I-l m— 1 

^ ] ^ ] Ik^m—i^i—k ^ ] ^ ] Il^iO/m—i—l ^ ^ ] Ii(X m —i 
k=0 i=k+l 1=0 i=l i=0 

:= Si + fi 2 + S 3 + B 4 + B B , m > 2. (4.8) 
Further we have 

m—2 m 

B 3 = ~E( E + E -^2)IkI m -iai-k 

k=0 i=k+2 i=k+l i=m 

m—2 m m—2 m—2 

— — E E Iklrn-itti-k ~ E + X] ^k^m-k 

k=0 i=k+2 k=0 k=0 

m-2 (m-2)-*; m -2 m-2 

= — X! E Ikha m -(k+l) ~ E Iklm-k-l&l + X! IkO-m-k, 
k=0 1=0 k=0 k=0 



m—2 m—l—2 

^4 = -E(E-E+ E )hha m -i-i 

1=0 i=0 i=0 i=m-l-l 

m-2 (m-2)-k m-2 m-2 

= — E E Iklia m -(k+l) + E IlQ-rn-l — E Illm-l-l&l- 
k=0 1=0 1=0 1=0 

Therefore the latter three terms in the right hand side of (|4.8|) becomes 

m -2 (m-2)-fc j 

-B 3 + B 4 + = — 2 E E hha m -(k+i) + - E hhF\. 

k =° 1=0 k+l=m-l 

k,l>0 

In this way, from Q4.8| ) we obtain 

^ m-l j m-2 (m-2)—*; j 

^m = — - E lilm-i ~ o E E hhF m -(k+l) ~ ~ E hhF\ 

1 i=l 8 fc=0 Z=0 , , , 1 

fc+f=m— 1 

fc,Z>0 

1 1 

= -oE J * J ™-* _ o E hIiF m -(k+i), m > 2, (4.9) 

2 i=i 8 

fc+Km-1 
A;,/>0 

by which we can determine any J m , m > 2, starting with Ji = — It is not 
difficult to find a homogeneous property among the terms of (|4.9| ). Thus we may 
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assume that 

m 

I m = J2 d n E F h ---F inl m>2. 



71=1 

i\ -\ \-% n =ra 



In general, the coefficients d n , 1 < n < m, should depend on m. But the following 
deduction implies that this assumption is possible. First from ( (4.9| ) we easily have 
^ 2 = ifs^i 2 ~~ I-^ 2 ' w hich leads to d 2 = t|o. When m > 3, ([Op becomes 



m—l m—l 
-^m = — - E Vm-i ~ g Fm ~ ~ E ^k F m-k ~ ~ E IkhF m -(k+l), (4-10) 
4=1 fc=1 fc+Z<m-l 

fc,Z>l 



in which the coefficients of the F™ yields the recursion relation (|4.6Q . In what follows 
we want to prove that I m , m > 0, determined above satisfy the relation ( |4. 1C| ) 



1 
J 

indeed. This may be shown by combining the following three equalities. First we 
have 

m— 1 m— 1 i m—i 

E hlm-i = EE^ E F h' ■■ F i k ^l d l E F h ' ' ' F ji 

1=1 7=1 fc=l . , , . . 1=1 . , , . 

nH h«fc=« JiH hj;=m-7 

«i,---,*fc>l Ji,— >i;>l 

m—l m—l m—i 

= EEEw E ^•••^ E I:, ■■■'■) 

k=i i=k i=i • , , • 

HH h7 fc =7 jiH Y3i=m-i 

ti,—,ife>i ii,— ,i;>i 

m— 1 m—k m—l 

= EEE« E ^•••^ E 

H H Nfc=« jiH hj;=m-« 

m 

= E E E ^•••■Pin- 

n=2 , , . . 

fe+i=Tl llH h7 n =m 

k,l>l ii,---,i n >l 

Similarly we can get the other two equalities 

m—l m 

^ ] Ik F m—k ^ ] d n —\ ^ ] ^ii ' ' ' Fi n , 

k=l n=2 

71 H h7n=77l 

7l>— i»n>l 
m—l 

E hhF m -(k+i) = E E E ^ii • • • -Pin+i ■ 

, ,^ „ 77 = 2 . . . . 

k+l<m— 1 t+]=n 7iH h7 n +i=m 

fe,Z>l M>1 7l,'",7 n -|-l>l 

Therefore the proof is finished. I 
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Theorem 4.2 If Pi and Qi, 1 < i < 3, solve the nonlinearized spatial system 
t/ien t/iere exist TV integrals of motion cti, < % < N — 1, of ( \3. 8j ) such that 



q = v / 2(< Pi, BQ 2 > + < P 2 , BQ 3 >), r = v / 2(< P 2 , 5Qi > + < P 3 , £Q 2 >) 
so/i>e i/ie following N-th order stationary AKNS equation 

N-l 

K N + a i K i = 0) 

i=0 



where < i < N, are defined by t\2.1Q) 



Proof: Noting that the expressions (|4.3| ) and (f4.4|) of b i+ i, c i+ i, % > 0, we can 
compute that 

N ~ ^ ( Ci x\ 



i=0 i=0 



h+i 



" * / v/2(< A*-iP 2 , SQi > + < A^'Pg, >) \ 
3 V v/2(< A i ^P 1 ,BQ 2 > + < A^P 2 ,BQ 3 >) J 



V <PXkW2k + <p2kty?,k 



i=0 j=0 k=l 

N N i 

fc=l i=0 jf=0 
AT )V )V 



01fc^2fc + 02fe^3fc 

02fc^lfc + 03fc^2fc 
<Plki>2k + </>2fc^3fc 



fc=l ;=o i=l 

Secondly we set 

N N N-l 



G(X) = IJ(A - Ai) =X) A A* = A w + 53 AA J 



i=l i=0 i=0 

Let us now choose 

N 

i=i 

which determines recursively 

N 

«at = Pn = 1, «z = A - 53 < / < JV - 1. 

i=Z+l 

due to Jo = 1. The a^, < i < iV — 1, are all integrals of motion of ( |3.8| ) since they 
are functions of Jj, < i < N . Further by G(\k) = 0, 1 < k < N, we see that 
Y,k=o a iKi — 0, which completes the proof. I 

The above theorem also implies that the potential determined by the Bargmann 
symmetry constraint (|3.5|) is a finite gap potential of the spectral problem (|2 . 1| ) . 
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Theorem 4.3 Under the control of the nonlinearized spatial system ( !Th ) , the non- 
linearized temporal systems (\3.Q) for n > can also be rewritten as the Hamiltonian 
systems 



Pit n — {Pi} Hn] 



B- 1 ^, Q Un = {Q l ,H n } = B- 1 ^, * = 1,2,3 (4.11) 



dQ l 



OR 



with the Hamiltonian functions 



H„ 



n 



4^0™ + ! 



E 



• • Fa 



ilH M m +i=rt+l 

«l,-",im + l>l 



where do = 1 and F m , m > 1, are defined by ( \3.1ty . 

Proof: We only prove the former equality of ( |4.11|) . We know that under the control 



of the nonlinearized spatial system ( |3.8| ), the results in Theorem ^4] holds. Hence 
we have 



P ltn = 2^^! + ^^^ 

i=0 i=0 

n i—1 n i—1 

= -2A n P x + 2^(]T J fc a,_ fe - I i )A n ~ i P 1 + v^E £ hk-kA^P* 

i=l k=0 i=l k=0 

n n—1 n n—1 n 

= -2 y £l k A n - k P 1 + y £ l h E ra^A^P. + Y.^ E V^k-kA^P, 

k=0 k=0 i=k+l k=0 i=k+l 

1 T p,-l dF l , 1 V- 1 r p-j djn-fc+1 1 A T r>-l dF n-k+l 

-LB + -^l k B — — — = ~ A l^ikB 



dQt 4 



k=0 



dQi 



k=0 



dQi ' 



where Oj, hi, i > 1, are given by (|3.18|). We further note the expression of I m , m > 0, 



defined by ( }4.5|) and then we may make the following performance 



P 



n—k+l 



1 r)F -ink 

^ IoB -'^±i + -yy / d m v F h ---F im B- 

4 dQi ^fi^ti . ^ 



ilH him=fc 

«l,-,'ra>l 



ilH Hm=k 



4 aQ x 4 ^ m + 1 



E 

ilH him+i="+l 

«l,-",»m+l>l 



P. ... P. 



1 8 n d 

4 



E F h - ■■ F im+1 = -B 



«iH Vi m +\=n+l 

ii,---,i m +l>l 



dQi 
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where we have accepted do — 1. The above manipulation is fulfilled for the case of 
n > 1. The case of n = 1 needs only a simple calculation. Thus the former equality 
of ( |4.il|) is true for n > 0. The latter equality of ([4.11 ) may be proved similarly 
The proof is completed. I 

The above theorem allows us to establish a sort of involutive solutions to AKNS 
soliton equations, which exhibits a kind of separation of variables for AKNS soliton 
equations. This is the following result. 

Theorem 4.4 The n-th AKNS soliton equation u tn = K n has the involutive solution 
with separated variables x, t n 

{ q = v/2(< <&<& n Pi(0, 0), Bg x H g^ n Q 2 (0, 0) > 

+ < <&<7%P 2 (0,0), £<£<& n Q 3 (0,0) >), 

(4.12) 

r = v/2(< <&<& n P 2 (0, 0), P^^Q^O, 0) > 

+ < ^ n P 8 (0,0), Bfl&flfe.Q 2 (0,0) >). 

where g y G denotes the Hamiltonian phase flow of G with a variable y and Pj(0, 0) 
and <5i(0, 0), 1 < i < 3, mm/ 6e arbitrary initial value vectors. 

Proof: Let 

P;(z, t n ) = g x H g*$ n Pi{0, 0), Q,(x, t n ) = g^Q^O, 0), 1 < % < 3. 

Then Pj(x,t n ) and Qi(x,t n ), 1 < z < 3, solve the nonlinearized spatial system 
( P-8| ) and the Hamiltonian system ( |4.11|) . However under the control of (|3.8| ), ( [4.11| ) 



is equivalent to the nonlinearized temporal system (|3.9| ). This shows that Pi(x,t n ) 
and Qi(x,t n ) also solve ( |3.8|) and ( p.9|) , simultaneously. Therefore the compatibility 
condition of (|3.8|) and (|3.9| ) is satisfied, i.e. ( [4. 12|) determines a solution to u tn = K n . 
In addition, since {H, H n } = 0, the Hamiltonian phase flows gfj, g^ may commute 
with each other. It follows that the resulting solution ( [4.12|) is involutive. The proof 
is finished. I 



5 Conclusions and remarks 

We have introduced a three-by-three matrix spectral problem for the usual AKNS 
soliton hierarchy and proposed the corresponding Bargmann symmetry constraint 
on this AKNS hierarchy. Moreover we have exhibited an explicit Poisson algebra 

{F 3 , l<j<N,F m ,m> 1} (5.1) 

on the symplectic manifold (R 6N ,lu 2 ) and further a binary nonlinearization pro- 
cedure is manipulated along with a sort of involutive solutions to AKNS soliton 
equations. When TV = 1,2, we have proved the nonlinearized spatial system ( |3.8|) 
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is integrable in the Liouville sense, indeed. But we don't know if we can take out 
enough independent integrals of motion among the Poisson algebra ( |5.1|) for a gen- 
eral integer N. We hope that this Poisson algebra suffices for proving complete 
integrability of the nonlinearized Lax systems. 

It should be pointed out that the Neumann symmetry constraint and the higher 
order symmetry constraints 

N r\ N r\ 

JT_i = J££rP, K m = JG m = JY,E j -^, (m>l), (5.2) 
j=i du j=i du 

may also be considered. These sorts of symmetry constraints are somewhat different 
from the Bargmann symmetry constraints because K_i is a constant vector and the 
conserved covariants G m , m > 1, involve the differential of the potential u with re- 
spect to the space variable x. In order to discuss them, we are required to introduce 
a new symplectic submanifold of the Euclidean spaces in the case of the Neumann 
constraint and new dependent variables, i.e. the so-called Jacobi-Ostrogradsky co- 
ordinates in the case of higher order constraints. Similarly, we can consider 
the corresponding r-symmetry (time first order dependent symmetry) constraints 
or more generally, time polynomial dependent symmetry constraints. Note that the 
similar Bargmann symmetry constraints have also been carefully analyzed for KP 
hierarchy |2(J and the symmetries in the right hand side of the Bargmann symme- 



try constraints are sometimes called additional symmetries |T0j and may be taken 



as sources of soliton equations [jn| . 



We remark that the finite dimensional Hamiltonian systems generated by nonlin- 
earization technique depend on the starting spectral problems. Therefore the same 
soliton equation may relate to different finite dimensional Hamiltonian systems once 
it possesses different Lax representations. AKNS soliton equations are exactly such 
examples. But we don't know if there exists an interrelation among the different 
finite dimensional Hamiltonian systems generated from the same soliton equation. 
In the binary nonlinearization procedure itself, there also exist some intriguing open 
problems. For example, why do the nonlinearized spatial system and the nonlin- 
earized temporal systems for n > under the control of the nonlinearized spatial 
system always possess Hamiltonian structures? We don't know either whether or 
not the nonlinearized temporal systems for n > are themselves integrable soliton 
equations without the control of the nonlinearized spatial system. These problems 
are worth studying in order to enrich integrable structures of soliton equations. 
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